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Abstract 

We apply the principles of discrete time mechanics discussed in earlier pa- 
pers to the first and second quantised Dirac equation. We use the Schwinger 
action principle to find the anticommutation relations of the Dirac field and 
of the particle creation operators in the theory. We find new solutions to the 
discrete time Dirac equation, referred to as oscillons on account of their ex- 
traordinary behaviour. Their principal characteristic is that they oscillate with 
a period twice that of the fundamental time interval T of our theory. Although 
these solutions can he associated with definite charge, linear momentum and 
spin, such objects should not he observable as particles in the continuous time 
limit. We find that for non-zero T they correspond to states with negative 
squared norm in Hilbert space. However they are an integral part of the dis- 
crete time Dirac field and should play a role in particle interactions analogous 
to the role of longitudinal photons in conventional quantum electrodynamics. 



1 Introduction 

Throughout this paper we shall use the acronyms CT to denote continuous time 
and DT to denote discrete time. The symbol T is used to denote our fundamental 
interval of time, which is assumed positive. The term CT limit and the symbol 
— s> will refer to the taking of T to zero and some integer n to infinity, such that 

nT — > t, where t is ordinary continuous co-ordinate time, assuming that this limit 

makes sense. This will frequently be the case, but it should always be kept in mind 
that DT mechanics is inherently richer in its range of possibilities than CT mechanics 
and includes trajectories for which such a limit is meaningless. This is particularly 
the case for the oscillon solutions which we shall discuss in the case of the DT Dirac 
equation. 
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This paper is the fourth in a series devoted to the construction of DT classical 
and quantum mechanics from first principles. This series of papers examines two 
questions: first, whether it is possible or consistent to rewrite all the laws of physics 
using discrete time rather than continuous time, and second, if such a programme 
were viable and self consistent, what would be the advantages (and perhaps disad- 
vantages) of such a discrete time formulation. Throughout this series we invert the 
standard view of discrete time. Rather than regarding it as a convenient approxi- 
mation to continuous time, suitable for say computer simulation, we suppose from 
the outset that our ordinary conception of time is itself an approximation, albeit an 
extraordinarily useful one, and that deep down there is a scale at which temporal 
phenomena occur in a discrete way and not in a differentiable or even continuous 
way. An analogy with water is useful here. For a vast range of phenomena, water 
may be sensibly modelled as a continuum, but sooner or later a molecular approach 
must be used. 

Our philosophy imposes a rigid discipline on our work. Our principles of discrete 
time mechanics must be as well specified as those of continuous time mechanics. Our 
laws of motion must be exact, and our invariants of the motion cannot be approxi- 
mately conserved but must be exactly conserved modulo those equations of motion. 
We should avoid introducing ad-hoc fixes for problems when they occur, but rely 
solely on the basic starting assumptions. In this respect we have found that our basic 
principles have guided us well, so that once a calculation has been started, we have 
found virtually no freedom of choice thereafter. What has emerged is a consistent 
theory of mechanics, both classical and quantum, which is genuinely different from 
continuous time mechanics, but which has sufficient overlap with it in so many ways 
as to make further investigation highly desirable. 

A question of some delicacy for DT mechanics is the issue of Lorentz covariance. 
Our DT mechanics is not manifestly Lorentz covariant. Fortunately, it turns out 
from our studies of scalar and Dirac field propagators that Lorentz covariance and 
the Poincare algebra appear to be broken in DT mechanics only at the order level 
so that our field theoretic scattering amplitudes are expected to have an expansion 
in powers of T such that the zeroth order terms are Lorentz symmetric. From our 
perspective, then, Lorentz covariance emerges as an approximate sjTiimetry of the 
complete mechanics, rather like isotopic symmetry is an approximate symmetry of 
nuclear physics. If as we imagine T is of the order of the Planck time or less, then there 
should be no problem in confronting current particle data, because special relativity 
has not been tested to anywhere near Planck scales. 

This raises an obvious question; why bother to discretise time if all that it does is to 
produce amplitudes which virtually duplicate conventional field theory? Our answer 
is that there are known problems with conventional formulations of field theory based 
on continuous time, such as divergences in Feynman diagrams, and our study is an 
exploration of an alternative approach to mechanics which might alleviate some of 
these problems. 

Our first paper, referred to as Principles I introduced basic principles for the 
temporal discretisation of CT classical and quantum particle mechanics. We have re- 
cently found a paper by Khorrami , who discussed various similar topics and which 
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is substantially in agreement with our results. Our second paper, referred to as Prin- 
ciples II 0, extended and applied our principles to classical field theory, including 
the Schrodinger equation, the Klein-Gordon equation. Maxwell's equations, gauge 
invariant electrodynamics and the classical Dirac field. The third paper. Principles 
III, [|] tackled scalar quantum field theory and there we discussed the construction of 
Feynman rules for DT scalar field theory. In that paper we found that otherwise hard 
or point-like vertices in CT ip^ Feynman diagrams are replaced by softened vertices in 
DT if^ theory, and that there is a conserved quantity analogous to energy in particle 
scattering processes, even though there is not a Hamiltonian in the theory. These 
papers should be consulted for further explanation of our notation, methodology, and 
motivation. In the present paper we apply the techniques of the previous papers to 
a discussion of the DT formulation of the quantised Dirac equation. We shall use 
natural units throughout, where c = h = 1. 

The next section is a review of the fundamentally important DT harmonic os- 
cillator, which provides the basic template for all our field theory propagators and 
generates the momentum space cutoff predicted in the free particle spectrum. Then 
we turn to the Dirac equation in one time and zero space dimensions. This provides 
a toy model which allows us to explore the novelties oi DT field theory. This model 
generates the normal spectrum of particle and anti-particle states, and also entirely 
novel solutions which, by virtue of their temporal behaviour, we call oscillons and 
anti-oscillons. The fundamentally bizarre property of oscillons is that they oscil- 
late in phase with a period twice that of the fundamental time T. We may readily 
appreciate that, in the temporal limit T ^ 0, these oscillon solutions cannot be 
accommodated within conventional CT mechanics. If T is extremely small but not 
zero, we may think of these oscillons as ghost like particles, capable of carrying mo- 
mentum, spin and charge, but not really a normal form of matter. It is an important 
result that our quantisation procedure, based on the DT Schwinger action principle, 
tells us that such single particle oscillon states, were they to be created, would be 
represented by state vectors having a negative inner product amongst themselves. 
These states would be regarded as unphysical, according to the usual interpretation 
of such vectors in state space. 

We then turn to the DT Dirac equation in one time and three spatial dimensions. 
The results are much as with the 1 + toy model discussed previously, but with 
some obvious differences. Now we can consider linear momentum and spin. We show 
how to construct the corresponding conserved variables using the Maeda-Noether 
theorem discussed in Principles I and //. We note that in DT we can examine field 
commutators and anti-commutators in fine detail, and there are one or two surprises 
with the Dirac field anti-commutators, related to the existence of oscillon solutions. 
We show also how to couple the DT Dirac equation to the DT Maxwell potentials 
in a gauge invariant way, and discuss the equations of motion and the conserved 
electric charge. Finally, we finish with a statement of the DT reduction formulae for 
the particle and anti-particle in and out states, in preparation for application to DT 
quantum electrodynamics, to be studied in subsequent papers of this series. 
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2 The Scalar field in 1+0 dimensions 



In this section we review the inhomogeneous DT harmonic oscillator equation, which 
serves as the basis for all DT free field equations, including the Klein-Gordon equa- 
tion, the Dirac equation, and Maxwell's equations. 

We showed in Principles I that the classical DT equation of motion for an har- 
monic oscillator in the presence of a source j is given by the Cadzow equation |^ 

PXn+l - 2aXn + pXn-1 = Tj^, P 0, (l) 

c 

where we use the symbol = to denote any equality holding on a true or classical 

c 

trajectory, and a and (3 are constants which are determined by the CT Lagrangian 
from which the DT system function is constructed. This equation can be written in 
the form 

{Un -2r] + U-^) Xn = J„, (2) 

c 

where r] = a/P, Jn = Tj^/P and Un is the classical step operator with the property 
Unfn = fn+1 for any variable / indexed by integer n. Equation (|) is a second 
order difference equation. We recall from our previous papers that for this particular 
equation there are two important dynamical regimes, called the elliptic and hyperbolic 
regimes, corresponding to r^^ < 1 and rj"^ > 1 respectively. These are separated by 
the parabolic barrier, rj^ = 1. Oscillatory behaviour occurs in the elliptic regime, 
which in field theories corresponds to the region where physical particle states occur. 
Singularities in propagators occur at the parabolic barrier m DT field theory, and 
avoiding these requires some careful discussion. 

To solve the equation of motion (|l]), first define the weighted differences 

^Xji = Xji ^jXji—Xi (3) 

and suppose that they satisfy the first-order difference equation 

Ax„+i = AAx„ + Jn (4) 

c 

for each value of n considered and for some complex-valued constants A and yU. Then 
we find 

A + = 2r/, A/x = 1, (5) 

which have solutions 



fx = r]± a/t^^ - 1, A = =F a/?7^ - 1. (6) 
In the case of forwards propagation, we find 

Xn — P ^ P + -P Jn 

c 

+... + PViv-2 + PV7v-i, N > n, (7) 

whereas for backwards propagation we find 

Xm = ~P Xn+l -\- P Xn -\- P Jn 

c 

+... + P'Jm+2 + P'Jm+i, M < n, (8) 
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where the polynomials P" = P'^{X,fi) are defined by 

P"(A, /i) = A" + A"- V + - + A/i"-^ + /x", (9) 
where (A, /i) = A + /i and P° (A, /i) = 1. From and we find 

pn-M-l^^ _|_ P^-'^-^xm = (^pn-M-lpN-n _|_ pN-n-lpn-M _ pi pn-M-1 pN-n-l-j^ 



n 

N~l 



_pn-M-lpN-n-lj^ ^ pn-M-l ^ pN-l-ij. (ig) 

i=n 

M+1 

p*-A/-l j.^ M < n < A^. (11) 



This is the desired relationship between the particle's position in the past (assumed 
known at time MT < nT), the particle's position in the future (assumed known at 
time NT > nT), and the particle's position at the present time nT. This equation 
incorporates Feynman (causal) propagation and the effect of the source terms and is 
an exact result. 

The next step is to look at the scattering limit oo, M — oo. Potential 

singularities at the parabolic barrier are avoided by taking the prescription t] 
7] + ie, e > 0, which corresponds to the Feynman m ^ m — is prescription in 
field theory, where e is infinitesimal and positive. It is convenient to introduce the 
parameter 9, related to r/ by r/ = cos^ and permitting it to take values in the complex 
9 plane only on the contour 9 = u — iv, where 

< u < 71, V = e > 0, (the elliptic regime) 

u = IT, 6 < V < OO (the hyperbolic regime) . (12) 



By writing rj = x + iy we see 



Vcoshw/ ~'~ Vsinhf 7 ^ 

Points in the elliptic regime lie on the upper half of an ellipse in the complex t] plane 
enclosing the points rj = ±1, which shows that the parabolic barrier is avoided. 
Turning to the P" functions, we use the fact that A/i = 1 to find 

PMA,/i) = ^^ /iVl. (14) 

jj, — jj, ^ 

Then taking /x = e*^ = e™"*"" we find 

lim — — - = —-^ : . (15) 



Now taking the scattering limit = — M ^ oo in (|TT]) we find 

oo 

-T Ay^j^, (16) 



m=—oo 
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where Xn^ satisfies the source free equation and the DT indexed Feynman propagator 
satisfies the equation 

(?7„ - + A^ = -r5„, r = r\ (17) 

where (5„ is the DT Kronecker deha defined in (|53D . Solving (0) we find 

Yp-i\n\e 

= 77—^, e = u-tv, (18) 
2i sm d 

which holds for the elliptic and hyperbolic regimes when we take into account the 
range of values the parameter 9 could take. 

Turning to particle theory and the Klein-Gordon equation, we saw in P] that the 
various constants in the above equations have the following parametrisations: 



where E = y p.p + m'^ is by definition the energy of the particle and p is its linear 
momentum. Now we define the transformed propagator 

oo 

A^(p,e)=T J2 e^"®A^(p), (20) 

n=— oo 

which satisfies the equation 

2(cos0 -r^) A^(p,0) = -rr. (21) 

Thence we find 

(p, 9) = -, (22) 

2 (cos (3 — 1] — le) 

using the Feynman +ie prescription for avoiding the singularities. In the above we 
have assumed G is real. We may define the propagator for complex values of B by 
analytic continuation of (^) . If now we introduce the variable po related to by the 
rule 

cos ® = — ^^9n (0) = sign (po) , (23) 

then we find 

1 T^o 

^ {Po ~ ^ + ie) ~^ 6 {pq — p2 — m'^ + «e) ' ^ ^ 

an exact result. From this we see the emergence of Lorentz symmetry as an ap- 
proximate symmetry of the mechanics. If Pq in the above is taken to represent the 
same thing in DT as it does in CT special relativity, i.e., the zeroth component of 
a four-vector, then the first term on the right hand side is clearly a Lorentz scalar, 
corresponding to the standard bosonic propagator of CT field theory. The second 
term on the right hand side is not Lorentz invariant but is down on the first term 
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by a factor proportional to T^. If as we expect T represents an extremely small time 
scale, such as the Planck time or less, then it is clear that the second term in (24) will 
be so much smaller in its effects than the first in general and so may be neglected in 
practice. This is why we expect our DT quantum field theory to be extremely well 
approximated by Lorentz covariant CT field theory in most situations. 

We can use the propagator and the DT Schwinger action principle to calculate 
the ground state expectation of DT time-ordered products of field operators, and 
hence extract the free field commutation and anticommutation relations. This was 
done in Paper III for the scalar field and we shall do the same for the Dirac equation 
in this paper. 



3 The quantised Dirac particle in 1+0 dimensions 

In this section we discuss the quantised DT Dirac particle in one time and zero 
spatial dimensions. This serves as the prototype for the Dirac field studied in the 
next section. We will use the Schwinger source function technique to obtain the 
ground state functional in the presence of the sources and from that we are able to 
extract the ground state expectation values of various anticommutators. 



3.1 Classical equations 



We start with the CT theory. Our dynamical variables are 'ip,^'^, which have two 
fermionic (anticommuting) degrees of freedom each: 



with = ip^'y^, where 

Then the CT Lagrangian is 

L = 



^1 

^2 



7 



1 
-1 



(25) 



(26) 



(27) 



where we assume the mass m is non-zero. Following Principles I and // we obtain the 
DT system function from the CT Lagrangian by considering the virtual paths 



A 



A, 



(28) 



and integrating TL{ip^, ip^) from A = to A = 1. With this prescription we find the 
time derivatives turn into differences. The result is the system function 



-- {2^+^.i7V„H-i + + ^;[+i7Vn + 2^;t7V„} , (29) 
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where k = mT. 

Next we add the Schwinger sources, which are taken as infinitesimal external 
disturbances to the system. This may be done in a number of ways. Because we will 
be interested in vacuum expectation values when the sources are switched off, it does 
not really matter what we choose to add and how. Our choice turns out to be the 
most convenient. We define the system function [rj] in the presence of external 
(i.e. non dynamical) fermionic sources 77, fj to be given by 

pn ^ ^ 1 y I -^^^ ^ r?„+lV^„+l + iJnVn + } (30) 

and then the Cadzow equation of motion 

^ =0 (31) 



gives 



or 



{^n+l - ^n-l) - y (^n+1 + 4V^n + ^n-l) = 'Vn. (32) 

^^i^n+i + 4k^„ + c^^„_i = 6rr/„, (33) 

c 

where ui = k + 3^7°, with a similar equation for the conjugate variable ip^ . 
Now define the non-singular matrix 



Cj = - — j- = e*^''' = cos ^ + i sin ^7 



\UJ\ 



e'^ 
e-^« 



(34) 



where \uj\ = ^/9~+~/? and the phase angle ^ satisfies (|^) . Then the equation of 
motion (RSp becomes 

6T 

{u-'Un - 27/ + uU-') = — r/„, (35) 

\uj\ 

where r] = —2k/\uj\. To solve this equation we first simplify the Dirac space depen- 
dence by the non-singular transformation = and then the equation of motion 
becomes 

6T 

{Un -2v + U~') 0„ = --u-'^Vn- (36) 

c \uj\ 

Using the results discussed in §2 we obtain the formal scattering solution 

oo 

0n = 0f -T J2 ^F-"'^-"'Vm, (37) 



m=—oo 



where ms a solution of the source free equation and the bosonic propagator satisfies 
equation (|I^ . This solution incorporates Feynman scattering boundary conditions. 
Transforming back to the original fields, we find 

oo 

^n=V^f Yl SF-'^Vm, (38) 



m=— oo 
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where 



on 
Op 



6 



2i\uj\ sin 9 



.^-i\n\e^n 



gg-i|n|e 



{cos + i sin 7°} 



2i\cu \ sin 9 

is the required DT Dirac propagator. It satisfies the equation 

6 



(39) 
(40) 



which is equivalent to 



Z7 



(C/„-C/-^)-^(C/„ + 4 + C/-^)U]^ 



(41) 



2T ^ " " ^ 6 
There are two important angles in this theory, given by 



(42) 



cos^ = 
cos,^ = 



-2k 

K 

+ 



sin 9 
, sin^ 



V9 - 3av2 

V9 + k2 
3 

V9 + k2 



(43) 
(44) 



where we assume we are in the elliptic regime 9 — 3k^ > 0. It is particular linear com- 
binations of these angles which characterise the solutions to the DT Dirac equation, 
as we shall now show. If we define the angles S = 9 — ^, a = 9 + ^ then we find 



cO _ 
Op — 



SI 



iV9-3K2' 
3 

iV9 - 3^2 
3 

iV9 - 3/^2 



e 



-inS 



e- 

gin^ 







-incr 



n > 



n< 0. 



(45) 



The significance of the angles 5 and ^ is the following. Assuming the fundamental 
parameter T is vanishingly small, a Taylor expansion about n — Q gives 



5 = 9-ic^K+—K^ + 0iK'), 
^ 180 ^ 

a = 9 + ^~n + ^K, + 0{K^) . 



(46) 



It turns out that the phase S is associated with physical particles in the CT limit 
T — > 0, whereas a is associated with wave-functions which change sign more-or-less 
once during each fundamental timestep T. This highly oscillatory behaviour is quite 
different to any behaviour encountered in CT mechanics and leads us to coin the 
term oscillon for solutions to the equations which behave in such a bizarre way. In 
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our theory there will be anti-oscillons as well as oscillons, in the same way there are 
anti-particles as well as particles. 

The reason oscillon solutions occur is not hard to understand. Our DT Dirac 
equation (|32D is, unlike the CT Dirac equation, manifestly second-order in nature. 
Such an equation will normally have two solutions. However, there are twice as many 
degrees of freedom when we take into account the Dirac space dimensions in the 
above theory, so we find a total of four solutions to the DT Dirac equation. These 
correspond to particle, anti-particle, oscillon, and anti-oscillon solutions. In Paper II 
we discussed the DT Schrodinger equation in exactly analogous terms, except we did 
not have a Dirac space structure. So in that example we found particle wavefunction 
and oscillon wavefunction solutions but no anti-particle or anti-oscillon solutions. 
We argued there that the oscillon wave-function solutions should not be accessible 
under ordinary conditions. This is supported by the result we discuss below, that 
quantum oscillon particle states have a negative norm and so are properly regarded 
as unphysical. 



3.2 Quantisation 

Having found the propagator the next step is to extract information from it about 
the field operators. We note that unlike CT mechanics, where we are usually able to 
impose canonical commutation or anticommutation relations directly, there is more of 
a problem in mechanics. This is because the concept of canonical conjugate mo- 
mentum becomes close to being redundant in our theory. Fortunately, the Schwinger 
source functional approach provides a powerful way of finding commutation or anti- 
commutation relations which are consistent with the quantum dynamics. 

Given a, DT action sum operator A'^^ [rj] , then for evolution from state M) 
at time TM to |a, A^) at time TN {N > M) in the presence of external sources, our 
DT Schwinger action principle is just 

5{a,N\p,M)r, = i{a,N\5A^^^ [rf\\j3,M)r^, N > M. (47) 

Here the action sum operator is given by 

N-l 

n=M 
N-l 

= ^{^'^ + k'^[Vn+li^n+l + Vni^n + '^n+lVn+l+'^nVn]}- (48) 
n=M 

{a,N\^J(3,M)^, N>n>M 
{a,N\^J(3,M)^, N>n>M (49) 
\{a,N\i;^\P,M)„ 



Then we find 



' -{a,N\P,M), 



i d 
—i d 



{a,N\P,M)r, 
{a,N\P,M), 
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^^(a,iV|/3,M), = i(a,iV|^^|/3,M)„ (50) 
^^(a,iV|/?,M), = |(«,iV|V;MlAM)„ 

^^(a,iV|/3,M), = i(a,iV|V^^,|/5,M),. (51) 

We may use these results to work out DT ground state time ordered products. For 
fermionic fields we note for example 

(0|T^;^^^|0) ^ {0\ijti^JO)Qn-m+'^{0\[ijt,k]mn-m (52) 

where T denotes the DT time ordering operator discussed in |l|, ||, ^, and 9„ and 5„ 
are the DT step and 6 functions respectively, with the properties 

6„ =1, n> 0, Sn = I, n = 0, 

= 0, n < 1, = 0, n ^ 0. (53) 

We note that in DT we can explicitly discuss what happens when the fields are 
at equal times, whereas in CT field theory, we tend to avoid this in discussions of 
time-ordered products, appealing to the properties of distributions. 
Now considering the quantum operator equation of motion 

-^-n 4k ^ 6T 

uj + —ip^ + uji>^_-y = —T]^, (54) 

\uj\ \uj\ 

we take its matrix element between the ground state, use the above relations (PU| — |^ 
and hence find the ground state functional 

Z[r^]^m, = exp\-tT' VrS^f'Vs] , (55) 



r,s=— oo 



suppressing Dirac space indices. Then using the result 

m^nai^mblO) = ^l^^^^^ ^^^^ 

explicitly showing the Dirac space indices, we find 

(0|T^„,^^,|0) = (57) 

or 

{Q\f^tJ^,\Q) = -2(5rV)6a. (58) 
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Using our results for the DT Dirac propagator in the previous subsection and choosing 
appropriate temporal indices m and n, we extract the following information: 



n>m: (0|V'+V'm|0) = 



V9 
3 


-3^2 


V9- 


3«;2 [ 
3 


3 


-3^2 


V9- 
3 


3/^2 [ 

e~ 


V9 - 
3 


3fi;2 


V9- 
3 


3«;2 [ 


V9- 
3 


3^2 


V9- 


3«;2 [ 


find 





-i(n-m)6'^rra-n^0 
_^-i{n-m)a g 



,Q A(n—m)B ,^,n—m 



^i{n—m)a 







^i{n—m)5 



i{m-n)e^m-n^O 

^i{n-m)S g 

g _^i{n-m)a 

^-i{n-m)S Q 

g _g-i(n-m)(7 



(0| i^n,i^n |0) 



37° 



2V9 - 3^2 

By using the operator equation of motion for the source free case 

4k 



\UJ\ 



we deduce 



(0K^J0) = -(0|^„^J0) 



37^ 



(59) 



(60) 



(61) 



(62) 



(63) 



(64) 



(65) 



V9 - 3fi:2 

We now use the above ground state expectation values and the reasonable assumption 
that free field anticommutators are c-numbers to write down the following anticom- 
mutators: 



{V'n+l>'0n} = 
{V'n>V'„+l} = 



6^7° ^ _i 
\uj\ 

-6i7° . 

, , CO, 



\UJ\ 



0, 



(66) 

(67) 

(68) 



knowing they are consistent with the operator equations of motion and with the 
propagator. This amounts to our quantisation prescription. 
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We note that the last result, (|68|) is surprising considering that the canonical CT 
equal-time anticommutator of Dirac field operators is non-zero, i.e. 



There are two points to be made here: 



Sab- 



(69) 



1. The Schwinger method in CT field theory does not lead to equal time anti- 
commutators directly. It deals with time-ordered products, which involve the 
Heaviside distribution (the step function) and then we have to consider temporal 
limits carefully. So for example the CT analogues of the results (|59| — 16^) are 



t>0: {0\^, {t)^,{0)\0) = 

{0\^t (0)^a (t)|0) 

t<0: (0|^,(0)V^^(t)|0) 

{o\ijAt)^t (o)|o) = 

from which we deduce 
lim 





e"^""* 



e^""* 

^imt Q 



g-fmt Q - 





ba 



ba 



ba 



(70) 



ba 



im {{0\^t (t) k (0) |0) + (0|^, (0) (t) |0) 



'ab. 



(71) 



which is consistent with 



2. The value zero in (^) is explained by the presence of oscillon solutions, which 
also manifest themselves in the relations (^ — . This suggests that the DT 
field operator 'ip^ is more complicated than its CT analogue ip (t) . 

3.3 The free particle Dirac equation 

We turn now to solutions to the free particle DT Dirac equation 



If we write ip^ = u)"^^ then 



(f/„ -2r] + f/-i) = 0, r] 



~2k 



(72) 



(73) 



Assuming this equation has solutions of the form 0„ = z^X) where x is a two compo- 
nent object and z is taken complex and non-zero, then we find 



{z-2r] + z-^) z^'x = 0, 



(74) 
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from which we deduce z = e^*^, where cos 6 = 7] satisfies . We will assume we are 
in the elliptic regime, so that 6 is real. Next, define the two component Dirac space 
basis functions 



(75) 





V2m 







u = 





V = 


\j2m 



Then solutions to ([7^ ) may be written in the form 

1 



2m 



(76) 



We note that there are four independent solutions, as discussed above. Hence finally 
we get 



iO (p. 

1 
2m 



f r —in& I * mtri , r i —ina , t* inSi 1 

|[ae +ce Jm + [cte +oe Jtij . 



(77) 



Second quantisation turns the coefficients a, 6, c, d, etc. into operators, so that the 
solutions to the quantised DT Dirac equation are 



1 

2m 
1 

2m 



ae 



m<5i 



{n+[a+e™^ + ce"*""] + i;+[d+e*"" + te-*"-^]} . (78) 



We may invert the relationship (assuming we are in the elliptic regime) to find particle- 
antiparticle creation and annihilation operators: 



2 sine 
—i 

2 sine 

i 

2 sine 
—i 



iO [ „inS 



T 



u^ip (t) e 
'if^ (t) ue 



imt 



-imt 



id I „,AnS 



2 sm e I 

and oscillon/anti-oscillon ladder operators: 



{t) ve 
► v~^^p (t) e' 



-imt 



(79) 



(80) 



2 Sine ^^n^^^'-^ne''^^^' 



2s-i^""^^-^^'^'-^"^ 



-i0 



d 



2 sine 
—i 



< 'p. 



f e 



2 sine I- ' ' " J ' T 

In the above equations we show what should happen in the limit 

— i> = lim{T — s> 0, n ^ oo, nQ t\ 



0. 



f82) 



(83) 
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assuming for example 

V'n+i-^W+^^W + 0(T2), etc (84) 

We note that we have not indexed the ladder operators with a temporal index n, but 
it is implied. These operators have the same temporal properties as the creation and 
annihilation operators found for the DT bosonic oscillator discussed in Principles I. 
Although such operators satisfy relations such as a„ = a„+i for example, they are 
not actually invariants of the motion in that they do not commute with the timestep 
operator In the terminology of Principles I, these operators are not compatible 
with the timestep operator. The resolution of this apparent paradox is to note that 
we are actually working in the Heisenberg picture. 

To understand further the basic properties of these operators, we consider the 
sixty four possible ground state expectation values of products of the eight operators 
a, a"*", Of these, only four are non-zero. We find 

(0|aa+|0) = (0|66+|0) 



V9-3k2' 

^ ^ — (S?77 

(0|cc+|0) = (0|drf+|0) = -=_. (85) 

V9 — 

This is consistent with the view that the operators a, b, c, and d annihilate the 
ground state ket |0) and that their adjoints annihilate the ground state bra (0|. 
If now we use the anticommutation relations (^ — ^8|) for the fields worked out 



previously, we find the following non-zero anticommutators: 

|a,a+) = \b,b+} = ^2m, (86) 

{c,c+} = |rf,d+| = ^^E^, ^-2m. (87) 

The fundamental result here is that the particle and anti-particle creation and anni- 
hilation operators lead to states with a positive inner product, whereas the oscillon 
and anti-oscillon operators lead to negative norm states, which means that oscillons 
and anti-oscillons should be regarded as unphysical. 

Finally, we find that the conserved charge Q operator, given by 

Q = Y^^t+i [3 - i^n + Y^i^t [3 + iKj°] v-^+i 

-^^n [3 - ^7°] - [3 + ^K7°] i^t (88) 



may be written in the form 

V9 - 3^2 



Q 



6m 



a^a — b^b + c^c — d^d 



^9) 



Taking the anticommutators into account, we see that particles, antiparticles, oscil- 
lons and anti-oscillons carry the same magnitude of charge. We expect therefore that 
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whilst oscillons and anti-oscillons would not appear as external in or out asymptotic 
particles, they can carry charge, linear momentum and spin, and should contribute 
in Feynman diagrams. In particular, there may be significant novel effects, with 
particles making virtual transitions to oscillons and back again. When we come to 
the field theory, discussed next, we expect to find that electric charge as a physical 
observable should only be associated with particles or antiparticles with momenta 
below the parabolic barrier. 



4 The Dirac equation in 1+3 dimensions 

Wc turn now to the Dirac equation in one time and three spatial dimensions. First 
we consider the free particle system. Using the virtual paths discussed previously 
the CT Dirac Lagrangian density 

£ (x) — ^i'li) {x) /dijj{x) — ^i-ii) (x) /dijj (x) — rmi) {x) ijj (x) (90) 
gives the system function density 

•^W = ^iW(x)Vn+lW-V':+l(x)V'nW} 



T I (x) H (x) + ^+ (x) H (x) 

6 \ +V;++i (x) (x) + 2V^+ (x) (x) 



(91) 



where 



^ =UV-cx-Ua-V +m(3 (92) 



is the standard Dirac one-particle Hamiltonian operator. The Cadzow equation for 
■0„ (x) is formally given by 



which gives the Cadzow equation for 



2T 



= {^„+i + 4^„ + ^„_i}, (94) 



where H = —icx- V + m(5. In this form we can readily identify the various parts of 
the non-covariant form of the CT Dirac equation. Another form is 

— AiT —>■ — 

A+V'n+l = ^^^V'n + A'V'n-l (95) 
c D 

where A+ = j + , A~ = ^ — . Likewise, we find the Cadzow equation for jp^ 

V.++,A- = —i^+H + V^+_iA+, (96) 



c 



6 
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where H = i V-a + mf3, A+ = ^ + , and A" = \- '^H . The equations in the 
form involving the A operators are useful for proving that various invariants of the 
motion are indeed constant in time. 
With the Fourier transforms 

^„ (p) = j ci^xe-P-^, (x) , V^:^ (p) = j ci=^xe^P->+ (x) , (97) 

we find 

A+ (p) i^n+i (P) = ^TH (p) ^„ (p) + A- (p) (p) , 

(p) A- (p) = -T^: (p) H (p) + V^ti (p) A+ (p) (98) 
c b 

where (p) = o: ■ p + A+ (p) = i + ^i/ (p) , and A^ (p) = i - ^^H (p) . Now 
consider the construction 

= / (p) (p) A+ (p) (p) + (p) A- (p) ^„ (p) } , (99) 



where C (p) is arbitrary. Using the equations of motion ( P^ ) we readily find 

C" = (100) 

c 

which shows that is a Logan invariant |^ of the system. 
4.1 The DT Dirac propagator 

In the presence of external fermionic sources the system function density becomes 

" ^ ^ " ^^^ + 2^1 (x) r^, (x) + (x) (x) / ' ^'"'^ 
which gives the Cadzow equation for ip^^: 

.fpn+l-'^n-l 

% ~ 



2T 

Likewise, we find the Cadzow equation for : 

AiT 
c 6 

Taking Fourier transforms and with the definitions 



i H {^„+i + 4V^„ + } - f3vn- (102) 



V^+^iA- = — V^+i/ + V^+_iA+ -zT r^+/3. (103) 



^^(p).^±|£|^, ..TE„i^(p).^ (104) 

the equation (|102|) becomes 

(p) f/„ - 2r/^ + (p) U-') (p) = TF^^ (p) Pf,^ (p) , (105) 
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where 



6 



-2k 



With our experience of the DT Dirac equation in 1 + dimensions, we may immedi- 
ately write down the formal solution 



(p) = i'f (p) - T 



r ™ (P) Vm (P) , 



(107) 



m=— CO 



where 'ip^^^ (p) is a solution to the free DT Dirac equation and the propagators 

5]^(p)=A^(p)r]"(p)^(p)7°, A^(p) = -— ( 
satisfy the equations 

n-\^)Un-2'nE + ^ij?)u-')~s'F{v) = -r£;i^ (p) 

(f/„-2r/^ + fZ-^)A^(p) = -VE5n. 



(108) 

(109) 
(110) 



Note that the bosonic propagator (p) in ( |108|) is dimensionless whereas the bosonic 
propagator in ( [18|) has the physical dimensions of a length in our system of units. 
To investigate the nature of the propagator, define the Fourier series transform 



~Sp (p, 6) ^ T ^ e^-''~Sl (p) 



:iiii 



where the parameter is taken real. Then we find 

n (cos e + 2) ^ - 3 sin el Sp (p, 6) = -3T/?. 



(112) 



We may now solve for Sp (p, ©) if we give a small imaginary term, according to 
the standard Feynman m ^ m ~ ie prescription. Hence we find 



-3T 



K (cos e + 2) - 3 sin e 



[/?2(cose + 2)2-9sin2e-ie] 



(113) 



In this form the propagator looks quite different to the standard CT propagator, but 
a suitable reparametrisation can change this. We introduce the parameter po (which 
should not be confused with Ep = a/p • p + m'^) related to the parameter B by 

6 - 2plT^ 



cos 



sign{Q) = sign (po) 



Then we find 
Sf ip,Q) 



pI — + ie p2 — + ie ^ ^ 



(114) 



(115) 
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From this we see that our propagator indeed looks exactly like the CT Feynman 
propagator for the Dirac field in lowest order in T. 



This result is important for two reasons. First, the reparametrisation (|1 14[ ) is 
precisely the same as the parametrisation (p3D used in the analogous expansion for the 
bosonic propagator, (P^. This means that the parameter fiowing through Feynman 
diagram networks has the same representation and interpretation for fermions as it 
has for bosons, and can be justifiably regarded as the DT analogue of energy, up 
to a factor of T. In Principles III we found that for scalar field theory, the sum of 
incoming G parameters was conserved at each vertex in a Feynman diagram, and 
so we conjecture that an analogous result holds for higher spin fields. The second 
important point is that this results shows us that Lorentz covariance for the Dirac 
equation emerges from DT mechanics at the same level of approximation as it does 
in bosonic theory. There is every reason, therefore to regard the loss of manifest 
Lorentz covariance in the theory as not a serious problem, provided that we work in 
the regime where T is close to zero in an appropriate sense. 



4.2 Field anticommutators 

With the Fourier series transforms implied by (|97D the action sum in the presence of 
external sources becomes 

N-l „ 
n=M 

N-l I, , 



Vn (P) 7°V^n (P) + (p) l^Vn (p) 



{Vn (p) 7°V^jv (p) + J (P) 7%iv (P 



{24 



+Vm (p) I'i^M (p) + i^M (p) i'vm (p)}- 

For the Schwinger action principle, we first define the functional derivatives 



6 



■Vm (y) = Sn-m5^ (x - y) 



6 



and 



so that 



5 



^Vn (P) 



6r]+ (x 
5 5 



(y) = Sn~mS^ (X - y) 



6ri+ (x) (p) 



(ixe 



^Vn (x) 



(116) 



:ii7) 



:ii8) 



5 



Sfln (P) 
6 



^+(q) = (27rf 5^ (p - q) 5. 



(119) 
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Then the Schwinger action principle gives 

^ -{a,N\P,M), = («,iV|7°V'„(p)|AM)„ N>n>M 



T Sfi^ (p) 

^^^-^(«,iV|/9,M), = {a,N\ijUp)l'\P,M)„ N>n>M (120) 
and so on. With this and Cadzow's equations of motion we find the vacuum functional 

= Z[0]exp j-.T^ ^ /7^^+(p)/^r"^(p)^^(p)|, (121) 



n,m=— oo 



where 



We may also write 



/ (^'^'"^^^ ^ ^""^ ■ ^^^^^ 



Z[rj]^Z [0] exp | -iT' E / ^^^^^ - y) } " (^^S) 

I, n,m=— oo J 

where the propagators 5"^ (x) satisfy the equations 

{^^^^^ + (^n + 4 + [/-^) } S-^ (x) = (x) . (124) 



Using the rule 



(P) ^mb (q) |0) = -i fe-" (P) /l ^ (27r)=^ 5=^ (p - q) (125) 

L J ba 

we find the following vacuum expectation values: 



WUla (P) ^nb 


(q)|o) = 


(9 + K2)sin^ 


(P)6a 


(27r)'53 


(P 


-q), 


(OlV'i (P) ^n+la 


(q)|o) = 


-18ie^^ ~_ 
(9 + K2)sin^ 


(P)a6 


(27r)'53 


(P 


-q), 


Wn+lb (P) ^na 


(q)|o) = 


18ie-^^ r 
(9 + K2)sin^ 


(P)6a 


(27r)'5' 


(P 


-q), 


(Ol^na (P)^n+16 


(q)|o) = 


-18ie^^ ~ + 
(9 + K2)sin^ 


(P)ab 


(27r)'53 


(P 


-q)- 



We notice that these are singular at the parabolic barrier. However, by taking anti- 
commutators and we arrive at the fundamental quantisation relations 

{V':-,i„(p),V'nJq)} = j^f^^k^ s^p - a) , 

{V'i(p),V'„+i„(q)} = (^|^A-(p)„,(27r)='53(p-q), (127) 

{V'i(p),V'„a(q)} = 0, 

which are, remarkably, free of any singularities at the parabolic barrier. If we had 
taken commutators instead, we would find that the singularities still occurred at the 
parabolic barrier. 
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4.3 Ladder operators 

Provided we are in the elliptic regime, the solution to the source free Dirac equation 
(q+ (p) V'n+i (p) - 2rjE^^ (p) + n (p) (p)) = (128) 

is given by 

1 2 

— ^ { [a (pr) e-^ + c (-pr) e^] u (pr) 



^n(P) 



+ [d+ (pr) e"^"^ + 6+ (-pr) e*"^] ^; (-pr)} . 



Then as outlined in the 1 + case, we find: 
for the particles: 



a{pr) = -^«+(pr){^„+,(p)e-^^-V'n(p)e^'} 



2 sine 

' 'K-,ii-p)e-'^-^U-p)e''}v{p 

(pr){V'n+i(-p)e^^-^n(-p)e-''} 



mi5 



b (pr) = 



2 sine 
— i 



r e 



in5 



—inS 



and for the oscillons: 



c(pr) 
c+ (pr) 
(i(pr) 
(pr) 



„in<7 



{^t+A-P)e'^-^U-P)e''}ui-pr) 
2 sin e ^"P"^ {^-+^ ^"P^ '"'^ " ^"P^ ''^'} 



2 sine 
— i 



-V'n (p)e-^^}^(-pr) 



2 sine 



Here we have used conventional Dirac momentum spinors defined by 



u 



(pr) 



Vm + E 



1 



Ur = 



(m- /») 
V (pr) = -u^ = 



1 



Vm + E Vm + E 



(m + E) X,. 

P '^Xr 
P <^^r 

(m + £;) 77, 



with 



(129) 



(130) 



(131) 



(132) 



(133) 



(134) 



(135) 



H (p) u (pr) = u (pr) , H (— p) (pr) = —v (pr) . 

Now using the field anti-commutators, we arrive at the following non-zero creation 
and annihilation relations: for the particles we have 



{a (pr) , a+ (qs) } ^ [b (pr) , 6+ (q.) } = ^y^^ ^rs {2nf 5' (p - q) 



(136) 
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whereas for the oscillons we find 



{c(pr) ,c+ (qs)} = {d{pr)J^ (qs)| = -^===Srs {27if 6' (p - q) . (137) 

All other anticommutators are zero. It is clear we should be in the elliptic regime 
for any of these anticommutators to make physical sense. Moreover, we see that 
even though their linear momenta may be in the elliptic regime TEp < \/3, oscil- 
lon and anti-oscillon particle states have a negative inner product and are therefore 
unphysical. This confirms the results of the previous section. 



4.4 Linear Momentum 



When there are continuous symmetries of the system function we may construct 
various invariants of the motion using the Maeda-Noether theorem discussed in Prin- 
ciples I [0, 0. First, consider the case when there is invariance under translation in 
space. From the results of Principles II 0] we construct the three components of the 
conserved linear momentum using the rule 



(ix 



d 



d 



d 



d 



138) 



This takes into account the anticommutation properties of the fields. Applying this 
rule to the translation invariant source free system function density ( pij ) we find 



(139) 



which has the appropriate limit 

lim^P" = j rfx ||#+ V^A - V^/'} (140) 
as expected. Using the equations of motion (^, ^6]) we readily find 

pn _ pn+1 (^4^^) 

c 

as expected. In the proof of this result we have to integrate by parts and assume 
that the wave-function falls off to zero at spatial infinity, which is to be expected for 
physical, i.e. normalisable states. 

Taking Fourier transforms, we find 



(ip 

(2^ 



p ^„ (p) A+ (p) (p) + (p) A- (p) ^„ (p) 



(142) 
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which on comparison with (|99|) shows that the hnear momentum is an example of a 
Logan invariant. After quantisation and using the resuhs ( 127 ) we find the commu- 
tators 



P",a+ (pr) 
P", c+ (pr) 



pc+ (pr) , 



P'^, 6+ (pr) 
P", d+ (pr) 



pb^ (pr) 
pd^ (pr) 



(143) 
(144) 



which shows that these operators do indeed create and annihilate excitations carrying 
definite linear momentum. 



4.5 Angular momentum 

Consider the free Dirac system function (pTI) . An infinitesimal rotation gives the 
following changes in the fields: 



where 



Sip^ = iu- J (5?/;+ = -iip^ J -uj, 



(145) 



(146) 



The system function is invariant to this transformation and so we may use the Maeda- 
Noether theorem to find the conserved angular momentum: 



With the Cadzow equations (^, ^6|) we find 

Jn ~ Jn— 1 

c 

as expected. We note in passing 

J =0. 



(147) 



(148) 



(149) 



a relation encountered in the CT Dirac equation. 



5 The charged Dirac equation 

We now discuss the coupling of electromagnetic potentials to the Dirac field. We 
follow here the DT formulation of the Maxwell fields discussed in Principles II 0, 
treating these fields as external, i.e., non-dynamical, and so in this paper we con- 
sider only the dynamics of the Dirac field. A full discussion of QED is reserved for 
subsequent papers in this series. 

DT gauge invariant electromagnetism, the electrostatic (scalar) potential 
is associated with the temporal links between successive instants of time, whereas 
the magnetic (vector) potential A is associated with the sites themselves. If 0„ (x) 
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denotes the scalar potential at position x on the temporal link {n, n + 1) and A„ (x) 
is associated with the position x at time n, then under a gauge transformation we 
have 



Xn+l (x) - Xn (x) 



T 



(f^'n (x) = 0„ (x) + 

AUx) = A„ (x) - Vx„ (x) , 

V'n(x) = X„(x)V'„(x), 



(150) 



where (x) = exp {—iqXn (x)} , q being the charge associated with the Dirac field. 

Turning to virtual paths, the presence of electromagnetic potentials requires a 
modification of the free field paths used previously. Following our discussion of the 
charged Klein-Gordon field in Principles II, we define the paths 



where Wn = exp {iq(f)^T) . The gauge covariant derivatives are given by 



(151) 



- iqAn^ ^„ 



T 



(152) 



where the virtual paths for the electromagnetic potentials are given by 

0„ = A„ = AA„+i + AA„. (153) 
Then under a gauge transformation, we find 



The gauge invariant system function for the Dirac field in external electromagnetic 
potentials may be given in the form 



F^=\T j dx {(V^+aV'n) + (V'KV'n)} 

where we use the angular brackets to denote integration over A, i.e. 

(/)^ f f{\)d\ 
Jo 

and 



(155) 



(156) 



On = iDn + iOt lDn " 
On = —iDn — iocl^n — 



(157) 
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Then we find the Cadzow equations of motion for the Dirac fields are 



(APF^^a^j + (Aiy-Aa-i^„-i) = 0, 



{^lOnXW-') + {^l_,On^,\Wti) = 0. (158 

c 



We note that, as found in Principles III, there is often an advantage in not eval- 
uating the integration over the virtual path until a late stage in a calculation. 
To go further, we introduce the following notation 

Al ^ \W-\^l, Bl^m^^t (159) 

Then we have 

^„ = v4„+i + 5„,, V'n=^^+l + 5n (160) 

and then the Cadzow equations give the relations 



(^n+l + 5„)) + (A^On-l {K + 5„_i)) = 



+ Bl) + + a-iA„) = 0. (161) 

These relations are very useful for proving charge conservation. Under gauge trans- 
formations where 

M/^ = X-iW^„X„, X„ = exp{-zgx„}, < = ^nV^„, (162) 

then 

A'„ = X„_iA„, 5; = X„5„, t^l = X„V'. (163) 

where 

= exp {-iq{\Xn+i + Ax J = X^^^x'l (164) 
To find the total charge consider the infinitesimal global gauge transformation 

^^n = -iQ^X^n^ ^i^t = k^X-iPt (165) 

and apply the Maeda-Noether theorem discussed in Principles II. Then we find the 
total charge can be written in the form 

Q" = iqT j dx - (E+O^A^+i)} (166) 

Then modulo the equations of motion ( |161| ) we readily find 

= g"-' (167) 

c 

and 



lim^Q" = q / rfx^/^+V (168) 



as expected. 
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5.1 Reduction formulae 

We anticipate now the construction of scattering amplitudes to be discussed in detail 
in subsequent papers of this series by giving the reduction formulae involving the 
Dirac particle and anti-particle creation and annihilation operators. If T denotes the 
discrete time ordering operator, discussed above and f/„ is the classical step function 
defined previously then we find the following reduced matrix elements: 



(ttoutloout (pr) (TO) |/3j„)i? 



K„,|(fO)a+ (pr)|/?J« 



\l^in)R 



\P^out 

\{TO)bl^ (pr) \(3in)R 



where 

In these matrix elements we have discarded the non-scattered components, as usual. 
In addition, we have not given any reduction formulae for the oscillon and anti-oscillon 
ladder operators. These we regard as not creating physically accessible states, and 
therefore, they should not be put on the same footing as the particle and anti-particle 
creation and annihilation operators. 

6 Concluding remarks 

The introduction of a non-zero T opens up a crack in OT quantum field theory which 
admits a number of new phenomena. In addition to the novelties encountered previ- 
ously in DT bosonic theories, such as a physical particle momentum cutoff, modified 
propagators and vertices, we find in the case of the Dirac equation the appearance 
of fermionic oscillon and anti-oscillon solutions. These should not survive in the CT 
limit on account of their extraordinary properties. For non-zero T they should cer- 
tainly participate in DT Feynman diagram processes as intermediate, virtual objects 
carrying charge, linear momentum and spin, but we expect them not to be observ- 
able as asymptotic in or out particles. We plan to investigate their role in DT QED 
scattering processes and regularisation in subsequent papers in this series. 



2 sin 6* 



J2 (P^) V{p){aout\f (p) 6) lAJ, 



n=—oo 

oo 



2 sine 



J2 (0^+ (p)) \PjV{^)u (pr) e--^ 



(169) 

oo 

^ E ('^out\f (-P) O) \(3jV{^)v (pr) 

n=— oo 

oo 

e-'^^'v^ (pr) V{-p){aout\f (o^n (-p)) \P^n) 

i=—oo 

(170) 



2 sine 



= n+(p)Un-2n+n{p)u-\ 



= Unn{p)-27i^ + u-'n+{p). 



(171) 
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